ЭЛЕКТРОДИНАМИКА ВЫСОКОЧАСТОТНОГО ФАКЕЛЬНОГО РАЗРЯДА 

И.А. Тихомиров, В.А. Власов, Ю.Ю. Луценко, А.А. Зорин 

Проведен анализ электродинамики высокочастотного (ВЧ) факельного разряда. Предложена и обоснована электродинамическая 
модель ВЧ факельного разряда. 


Введение 

ВЧ разряды, в соответствии с установившейся 
классификацией [1-4] называются ВЧ-индукцион- 
ными (ВЧИ), ВЧ-емкостными (ВЧЕ) и ВЧ-дуговы- 
ми (ВЧД) разрядами. ВЧ-индукционный, так на- 
зываемый безэлектродный Н-разряд, был описан 
еще в 1884 году [5]. В ВЧ -дуговых разрядах между 


электродами внутри разрядной камеры возбужда- 
ется дуговой Е-разряд. Обычный ВЧЕ разряд воз- 
буждается в разрядной камере при наложении 
электродов на внешние диэлектрические стенки 
камеры. Среди ВЧ разрядов особое место занимает, 
так называемый, ВЧ факельный разряд. Разряд был 
открыт и так назван С. И. Зилитинкевичем в 1928 г. 
[6], т.к. разряд (рис. 1) с электрода - 1 горит в виде 
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факела, в котором высокотемпературный канал - 2 
окружен светящейся диффузной оболочкой - 3, ди- 
аметр которой значительно больше, чем у канала. 
Энергия к разряду подводится от ВЧ генератора че- 
рез электрод - 1. 

По своей природе и свойствам ВЧ факельный 
разряд является одноэлектродным Е-разрядом, ко- 
торый горит с учетом емкостной связи факел-зем- 
ля [7]. Вид, форма, режимы горения факельного 
разряда, его специфика описаны в работах [1-3, 8- 
19]. Разряд возбуждается в любой (включая и водо- 
род) газовой среде при давлениях 10...5-10 5 Па. Ис- 
следования показали, что плазма ВЧ факельного 
разряда обладает ярко выраженной неравновес- 
ностью. В зависимости от условий электронная 
температура в плазме разряда изменяется от 5000 
до 35000 К и выше, а газовая (молекулярная) - от 
1500 до 5000. ..9000 К. В каждом конкретном случае 
электронная температура в несколько раз (иногда 
на порядок), выше, чем газовая. Степень иониза- 
ции плазмы может изменяться с учетом электрон- 
ной плотности ІО 15 ... ІО 22 м 3 в диапазоне значений 
10 ' 8 . .. 10 1 . 

В свое время [18, 20] была предложена прибли- 
женная электродинамическая модель разряда. Ка- 
нал разряда рассматривается в этой модели как од- 
нородный проводник с распространяющейся по 
нему основной ТМ-волной. Выделяемая в разряде 
энергия определяется омическим сопротивлением, 
а ее рассеяние - радиальной теплопроводностью. 
Моделирование для более мощных ВЧ факельных 
разрядов, когда их геометрические размеры стано- 
вятся соизмеримыми с длиной двухпроводной ли- 
нии анодного контура ВЧ генератора, рассмотрено 
в работах [1, 21, 22]. С учетом слабого скин-эффек- 
та диссипация электромагнитных волн различного 
типа в разряде, а не только основной ТМ волны 
проанализирована в [ 17]. В плазме канала, со слабо 
изменяющейся в осевом направлении величиной 



Рис. 1. ВЧ факельный разряд 


комплексной диэлектрической проницаемости, 
могут возбуждаться симметричные поперечно- 
электрические (ТЕ) и поперечно-магнитные (ТМ) 
волны, а также несимметричные волны, являющи- 
еся гибридными (рис. 2). Характер электромагнит- 
ного поля разряда изучался (рис. 1) с помощью од- 
новитковой рамки - 4, присоединенной к осцил- 
лографу. Численные расчеты и эксперименты [17] 
показали, что при средних мощностях определяю- 
щий вклад в энергобаланс разряда вносит основная 
ТМ волна. Вклад побочных ТМ волн незначителен 
и слабо зависит от мощности и частоты. Следует 
учитывать, что вклад ТЕ волн при средних мощ- 
ностях (до 10 кВт) меньше, чем от основной ТМ 
волны, но с ростом мощности разряда и частоты 
доля вклада от ТЕ волн возрастает [17]. 

В работе [23] изучено распределение (рис. 1 и 
рис. 2) компонент электромагнитного поля Н ' и Е г 
и их затухание в канале разряда; определено ради- 
альное распределение удельной электропровод- 
ности плазмы в разряде. Перегиб профиля удель- 
ной электропроводности наблюдается на границе 
раздела "канал-диффузная оболочка". При этом 
максимум магнитного поля Н наблюдается на пе- 
регибе, а максимум Е г - в пределах диффузной обо- 
лочки разряда - 3. 

Данные о свойствах и особенностях ВЧ факель- 
ного разряда, специфике его электромагнитного 
поля, логично позволяют нам перейти к анализу 
электродинамики разряда с обоснованием новой 
электродинамической модели разряда в отличии от 
модели, предложенной в работах [18, 20]. 

1. Распределение электрических токов 

в канале высокочастотного факельного разряда 

В отличие от работ [18, 20] будем рассматри- 
вать канал разряда в виде излучателя конечной 
длины. В этом случае канал разряда можно предс- 
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Рис. 2. Распределение компонент электромагнитного поля в 
канале разряда 


22 



Естественные науки 


тавить [25] в виде тонкого линейного излучателя 
конечной длины. Распределение токов в таком из- 
лучателе в соответствии с интегралом Кирхгофа- 
Гюйгенса [25] будет зависеть от Ну. 
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Здесь: /(/) - распределение электрических то- 
ков вдоль линейного излучателя; / - расстояние от 
элемента излучателя с координатой і до точки 
наблюдения Р, имеющей координаты х, у, % к - 
постоянная распространения электромагнитной 
волны в непроводящей среде; Н - аксиальная ком- 
понента магнитного поля (рис. 1). 

Электрическую составляющую электромагнит- 
ного поля можно легко определить через извест- 
ную величину магнитной составляющей посред- 
ством следующего выражения [25]: 


Е = ( 1/ Дѵе)ѴхЯ. 


Прежде, чем проводить анализ эксперимен- 
тальных результатов, рассмотрим распределение 
компонент поля для случая "электродинамической 
модели" разряда по работе [18]. Согласно этой 
"электродинамической модели" вдоль канала 
ВЧФР распространяется затухающая поперечно- 
магнитная волна. При этом предполагается, что от- 
ражение ее на конце канала разряда несуществен- 
но, и им можно пренебречь. В этом случае распре- 
деление электрических токов в канале ВЧФР запи- 
шется в виде [18]: 

і= ѴД у = а +у А (2) 

где у - комплексный коэффициент распростране- 
ния электромагнитной волны (волновое число) 
вдоль канала разряда. 

Преобразуем интеграл (1) с учетом выражения 

( 2 ). 
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где: 11= е аі 'со&[-(Рі' + кг 1 )]; 

Ѵ= е‘ аг ’8Іп[-(/Д' + к/)]. 

Результаты расчетов величин интеграла (3) и, 
соответственно, распределения аксиальной компо- 
ненты магнитного поля вдоль оси ВЧФР для раз- 


личных значений у 



при у =3- 10 2 м приведены на рис. 3. Расчет прово- 


дился методом Симпсона. 

Как видно из рисунка, постоянство аксиальной 
компоненты магнитного поля для случая "электро- 
динамической модели" [18] может наблюдаться 
только на расстояниях от электрода, не превышаю- 
щих нескольких сантиметров, что составляет вели- 
чину не более 15. ..20 % от длины канала ВЧФР. За- 
метим также, что и характер изменения фазового 
сдвига вдоль оси ВЧФР, рассчитанный на основе 
"электродинамической модели” не соответствует 
полученной нами экспериментальной кривой. 

Таким образом, полученные нами эксперимен- 
тальные результаты не могут быть интерпретирова- 
ны соответствующим образом в рамках "электро- 
динамической модели" [18]. 

Перейдем теперь непосредственно к задаче оп- 
ределения распределения электрических токов в 
ВЧФР по результатам измерений аксиальной ком- 
поненты магнитного поля в ближней зоне его излу- 
чения. 

В общем случае задача нахождения распределе- 
ния электрических токов в излучателе по данным 
измерений создаваемого им поля достаточно слож- 
на, и ее решение математически разработано [25] 
только для анализа компонент поля в дальней зоне 
излучения. Однако ее приближенное решение нес- 
ложно получить, основываясь на том или ином 
численном методе расчета интеграла Кирхгофа- 
Гюйгенса. 

В соответствии с методом трапеций прибли- 
женное значение интеграла (2) запишется в виде 
[26]: 
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Подынтегральную функцию Д/д';/) можно 
представить в виде: Д /;/;/) =а(гд')|у= СОП5( Ч(/), 
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Выбрав соответствующим образом п экспери- 
ментально полученных значений аксиальной ком- 
поненты магнитного поля мы можем записать п 
линейных уравнений на п неизвестных величин. В 
матричном виде данная система уравнений запи- 
шется следующим образом: 

И ■ И] = [б], (4) 

где [ Д - вектор-столбец неизвестных электричес- 
ких токов, соответствующих определенным участ- 
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кам канала разряда; [А] - матрица коэффициентов 
а(^х'); [ В\ - вектор-столбец экспериментально по- 
лученных значений аксиальной компоненты маг- 
нитного поля. 

В качестве примера рассмотрим решение зада- 
чи о нахождении распределения электрических то- 
ков излучателя по известным значениям аксиаль- 
ной компоненты магнитного поля для случая 
ВЧФР с длиной канала 1=18 см, что соответствует 
мощности, вкладываемой в разряд \Ѵ «2 кВт. 

Разобьем канал разряда на 6 отрезков длиной 
3 см. В этом случае матрица коэффициентов а(х',І) 
при у =30 мм будет иметь следующий вид: 
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Здесь с целью упрощения мы пренебрегаем за- 
паздыванием сигнала, что выражается в отбрасы- 
вании мнимых частей коэффициентов, величина 
которых для ближней зоны излучения сравнитель- 
но мала. 

Вектор-столбец значений аксиальной компо- 
ненты магнитного поля составим, пользуясь имею- 
щимися экспериментальными данными. В данном 
случае он будет иметь следующий вид: 

1,00 + у 0,00 

0,99 + у 0,09 
В = 

0,90 + у 0,19 

0,96 + у 0,26 


"Выбросы" амплитуды электрического тока в 
крайних точках связаны с краевыми эффектами, а 
также с приближенностью метода решения задачи. 
Однако в реальности краевые эффекты не столь су- 
щественны для случая факельного разряда, так как 
с одной стороны канал разряда электрически свя- 
зан и является своего рода продолжением электро- 
да, с другой стороны канал разряда не имеет резкой 
границы. Расчет, проведенный с учетом диффуз- 
ной части разряда, подтверждает вышеприведен- 
ные соображения. 

На рис. 4 приведены распределения электри- 
ческих токов вдоль оси свободного ВЧФР с учетом 
и без учета диффузной части разряда, а также для 
случая экранированного ВЧФР. 

Как видно из рис. 4, постоянству амплитуды ак- 
сиальной компоненты магнитного поля в ближней 
зоне излучения разряда соответствует распределе- 
ние электрического тока с неэкспоненциальным 
характером изменения амплитуды вдоль оси разря- 
да. 

Заметим, что отсутствие существенных измене- 
ний электрических параметров плазмы разряда 
вдоль его оси можно объяснить наличием: 

- слабого затухания электромагнитной волны, за- 
питывающей разряд; 

- отраженной электромагнитной волны в конце 
разряда. 

При отсутствии затухания электромагнитной 
волны горение разряда невозможно, так как в этом 
случае не происходит передачи энергии электро- 
магнитной волны к плазме разряда. 

Такими образом, с целью интерпретации полу- 
ченных расчетных и экспериментальных результа- 
тов мы должны предположить наличие отражен- 
ной электромагнитной волны в канале факельного 
разряда [27-29]. 


0,94 + у 0,32 
0,82 + у 0,57 
0,69 + у 0,72 


Здесь с целью упрощения анализа мы использу- 
ем относительные, а не абсолютные значения ак- 
сиальной компоненты магнитного поля. Переход к 
абсолютным значениям несложно осуществить ум- 
ножением полученных значений амплитуд токов 
на соответствующий множитель. 

Решение системы уравнений (4) найдем путем 
использования метода Гаусса. После преобразова- 
ний получим значения величин электрических то- 
ков, соответствующих определенным участкам ка- 
нала ВЧФР. Вектор-столбец электрических токов 
будет иметь вид: 
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Рис. 3. Расчетное распределение аксиальной компоненты 
магнитного поля вдоль оси ВЧФР: 1) относительная 
амплитуда ; 2) фазовый сдвиг. (1, 2) - /у/ = 45 м ’; 
атду =7 1/4 ( Т, 2 ) - ІуІ=30 м’; агду =п /4 
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Рис. 4. Распределение амплитуды (1) и фазового сдвига (2) 
высокочастотного тока в канале ВЧФР: без учета (а) и 
с учетом диффузной части разряда (б); (в) - 
экранированный разряд 


2. Электродинамическая модель высокочастотного 
факельного разряда с учетом отраженной 
электромагнитной волны 


Впервые высокочастотный факельный разряд 
был рассмотрен в виде линии с распределенными 
параметрами конечной длины в работе [30]. Авто- 
ры вышеуказанной работы справедливо указали на 
недостатки электродинамической модели ВЧФР 
[18], основанной на представлении разряда в виде 
цилиндра бесконечной длины, которые выража- 
лись, прежде всего, в расхождении теоретических 
расчетов с экспериментальными результатами. За- 
метим, что модель бесконечного цилиндра приме- 
нима только в случае отсутствия в канале ВЧ фа- 
кельного разряда отраженной электромагнитной 
волны. Другими словами, учет конечности длины 
канала ВЧФР при теоретическом рассмотрении 
выражается прежде всего в учете наличия в канале 
разряда отраженной электромагнитной волны, что 
и сделано в предлагаемой модели. 

Рассмотрим волну электрического тока, расп- 
ространяющуюся вдоль линии конечной длины 
(см. рис. 5). В общем виде ток, протекающий в точ- 
ке і, определяется [26] следующим образом: 




У2 ]С, ~У 2 
е‘ + се е 


(5) 


Здесь: /, / 0 - комплексные амплитуды ВЧ-то- 
ков, соответственно в точках і и і = 0: у = а + //3 - 
комплексный коэффициент распространения 
(волновое число) электромагнитной волны в кана- 
ле ВЧФР; се & - комплексный коэффициент отра- 
жения электромагнитной волны на конце линии. 

Выясним далее, при каких значениях парамет- 
ров а, Р, С, с, рассчитанное на основе эксперимен- 
тальных результатов распределение электрическо- 
го тока в канале ВЧФР может быть удовлетвори- 
тельно описано выражением (5). 

Результаты расчетов распределения амплитуды 
и фазового сдвига электрического тока вдоль кана- 
ла ВЧФР, проведенные согласно выражению (5) 
для различных значений а, Р представлены на рис. 
6. Сопоставляя кривые распределений электричес- 
кого тока, изображенного на рис. 4, можно видеть, 
что наибольшее совпадение распределений элект- 
рического тока, полученных, соответственно, на 
основе экспериментальных данных и в результате 
расчетов по формуле (5) наблюдается в случае 
а = 3 м'\ Р< 30 м~\ % = 0, с = 1. 


2 

I 4 
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Рис. 5. Схема линии конечной длины 
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Малую величину коэффициента затухания 
электромагнитной волны, распространяющейся 
вдоль канала ВЧФР, подтверждают также измере- 
ния фазового сдвига между радиальной компонен- 
той магнитного поля Е г и аксиальной компонентой 
магнитного поля Н г Заметим, что для цилиндри- 
ческих структур в непроводящей зоне отношение 
комплексных амплитуд Е г и Н ' непосредственно 
определяется [3 1 ] величиной волнового числа к. 


Е г Іі/лсо 



( 6 ) 


Учитывая, что величина - действительная, 

к I 

из выражения (6) получим: 




/агсТ ё 

Се 


1 

а 


(7) 


где С - действительная константа; /3 - коэффициент 
затухания; а - фазовый коэффициент (И = а +//?). 

Таким образом, согласно выражению (7) по ве- 
личине фазового сдвига между радиальной компо- 
нентой электрического поля и аксиальной компо- 
нентой магнитного поля однозначно определяется 
отношение коэффициента затухания к коэффици- 
енту фазы электромагнитной волны, распростра- 
няющейся в плазме ВЧФР. 

Эксперимент по определению характера маг- 
нитного поля ВЧФР позволяет также провести 
оценку разности фаз между компонентами поля Е г 

ассмотрим более подробно смещение кривой 
сигнала I в зависимости от угла вращения рамки <р 
(рис. 1), обусловленное наводками компоненты 
поля Е г Сигнал от рамки, поступающий на вход 
измерительного прибора определяется как сумма 
ЭДС, наводимых магнитным и электрическим по- 
лями 



р/а 



Рис. 7. Зависимость соотношения [3/а от мощности разря- 
да И/ 


1 = 8 , 


+ е 


Н, 
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Учитывая, что геометрические размеры рамки и 
частота электромагнитного поля постоянны, мож- 
но записать: 


Іе 




= Ае 


/Ч / р № н 

Е '- + Ве Ч> 


81П ф. 


где А и В - постоянные действительные величины, 
характеризующие величины сигналов, поступаю- 
щих на вход измерительного прибора соответ- 
ственно от электрической и магнитной компонент 
поля. 

Преобразуя выражение (7), получим: 


А 2 +В 2 зіп 2 (р + 


_Л 

Е г 


+2АВ сов 


-4 х 


н, 


(р 


8Ш (р 


( 8 ) 


ЭзтЧ' 27 + І?81пЧ ' Е7 8 ПК» 

Ч 7 , = агс!§ ті_ (9) 

1 ЧсозЧ* г + ВсокѴ „ 8 Іпю 

Ь г Пу 

Пусть для определенности Ч 1 = 0 • 

І-лу 

Тогда выражение (9) примет следующий вид: 


ЯзіпЧ' 


А* і = агс!§ 


тт 8Ш (О 


А + В С08 Тц 8ІП (О 


При (р = 0 очевидно, что Ч^ = 0. Получим: 

С 


с!§4 / 7 = 


8ІП Ч* 


тт 81П (р 

Л(р 


(10) 


Здесь С, = — . 

О 

п ПК 

Выражение (10) в случае (р = — и (р = 

2 2 

запишется в виде системы уравнений: 

= _с 

2 


с!§Ч / і 


т _п ~ . ф + с ^ н > С 11 ) 

8Ш Ч* „ п (р 


н, 


<р 


сі§Ч / і 


<р 


= _7 Т - 


С 


(р 8ІП Ч^ 


•( 12 ) 


(Р 


Для случая ВЧФР, горящего при атмосферном 
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давлении мощностью \Ѵ=\ кВт на расстоянии от 
электрода / =60 мм, имеем: 




I 



Ч> 


(р 


-_п_ 


= -40 


Соответственно, в результате решения системы 
уравнений (11-12) получим, ^ тт = 16 ; 

п (р 

С-122. р 

Используя (7), будем иметь — = 0,287. 

а 


Результаты измерений 


I 

а 


в 


зависимости от 


мощности разряда приведены на рис. 7. 

Как видно из рис. 7, величина — находится в 

а 

пределах 0,25. ..0,35. 

Учитывая, что длина электромагнитной волны 
имеет величину большую, чем длина канала, т.е. 
2л/ а > X, получим, что падение интенсивности 
поля в е раз происходит на расстоянии от электро- 
да, которое сопоставимо или превышает длину ка- 
нала разряда, так как / 3 < 1,8/Х. Таким образом, по- 
лученная нами оценка величины коэффициента 
затухания для ВЧФР мощностью 0,7... 1 кВт в нес- 
колько раз меньше соответствующей величины, 
приведенной в работах [18, 20], что говорит в поль- 
зу предлагаемой новой модели разряда. 

Рассмотрим распределение мощности диссипи- 
рующей в плазме ВЧФР электромагнитной энер- 
гии. 

Мощность джоулевых потерь, возникающих в 
результате диссипации электромагнитной волны, 
запишется следующим образом [20]: 

УѴ = 2л-} } К |/|” гсігсіг . (13) 

оо 


Здесь |/| - действующее значение высокочастот- 
ного тока; г - радиус канала ВЧФР; Я - величина ак- 
тивной составляющей сопротивления плазмы 



0,25 0,5 0,75 1,0 

2=2/Ь 

Рис. 8. Зависимость сІѴѴ/сіг от длины канала разряда 


ВЧФР, приходящейся на единицу длины, которая в 
случае слабого скин-эффекта определится следую- 
щим образом: 


К = (рлг 2 ) Л . (14) 


Используя выражение (5) для высокочастотно- 
го тока, распространяющегося по линии конечной 
длины, получим выражение для действующего зна- 
чения тока: 


/ = /„ 


е ( ' хг соз (5і + 

+Се С08 (ср-^г) 


-і2 


+ 


+ 


е ал $іп I Зг + 

+Се аг зіп(<р -/Зг) 


• (15) 


С учетом выражений (14) и (15) интеграл (13) 
для случая С = 1, (р = 0 запишется следующим об- 
разом: 

УѴ = 2сг _1 1/ 0 |~ } (скіаг -зіп 2 (Зг)сІ 2 . (16) 
о 

Мощность, диссипирующую в канале ВЧФР, 
найдем путем интегрирования выражения (16). 

Получим: 

ЩЬ) = 


= 2(т 1 



сН2аЬ 

+ 

2 а 

+ — — зіп ВЬ соз ВЬ — — 
4 р 2 


(17) 


Найдем распределение плотности источников 
диссипации электромагнитной энергии вдоль ка- 
нала ВЧФР Мощность, приходящаяся на единицу 
длины канала ВЧФР, характеризуется величиной 

Ш 

. Из выражения (17) следует, что: 

СІ2 

сПѴ _ _1 . ,2 , , _ . 2 \ 

= 2 сг / 0 (ск2а.2 -зпГ (Зг). 

СІ2 

р ^УѴ 

Результаты расчетов величины для случая 

СІ2 


разряда с длиной канала X приведены на рис. 8. 
Как видно из рис. 8, по длине разряда можно выде- 
лить три зоны интенсивного энерговыделения - на 
краях и в центральной части канала разряда. 

Подводя итог вышесказанному, делаем такие 
заключения: 

1. Полученное на основе теоретических расчетов 
и экспериментальных результатов распределе- 
ние электрического тока вдоль канала ВЧФР в 
отличии от работ [18, 20] носит неэкспоненци- 
альный характер. 
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2. Расчеты, проведенные на основе предложенной 
модели разряда, представленного в виде линии 
конечной длины с распределенными парамет- 
рами, подтверждают предположение о наличии 
в канале ВЧФР отраженной электромагнитной 
волны. 

3. Отражение электромагнитной волны и тока в 
конце канала разряда носит характер, близкий к 
синфазному. 

4. Длина электромагнитной волны, распространя- 
ющейся вдоль канала ВЧФР, как правило, пре- 
вышает длину канала разряда, т.е. выполняется 
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